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^- , Abstract. Let V be a vector field distribution on manifold M . We give an efficient algorithm for 

^"^ ' the construction of local coordinates on M such that V may be locally expressed as some partial 
prolongation of the contact distribution Cq , on the 1*^* order jet bundle of maps from K to K', 

\^ ' g > 1. It is proven that if V is locally equivalent to a partial prolongation of Cq then the explicit 

^^ , construction of contact coordinates algorithmically depends upon the determination of certain 

• ' first integrals in a sequence of geometrically defined and algorithmically determined integrable 

i -^^ I Pfaffian systems on M . The number of these first integrals that must be computed satisfies a 

J^ , natural minimality criterion. These results therefore provide a full and constructive generalisation 

of the classical Goursat normal form from the theory of exterior differential systems. 
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O ■ 1. Introduction. 

o 



Recent years have witnessed progress in the creation of geometric characterisa- 
tions of contact distributions over jet bundles [1,9,4]. The problem of obtaining 
a simple, invariant characterisation of the corresponding partial prolongations is a 
natural one and has apparently remained open, even in the case of partial prolon- 
gations of jet bundles of maps from the line. That is, partial prolongations of the 
k>( ; contact distribution Q ' over the first order jet bundle of maps M ^ M'^', g > 1. 

i_( ' In [6], a solution to the characterisation problem for partial prolongations of 

- - ■ Cq was presented. Specifically, simple geometric conditions, expressed in terms of 

the derived type of V C TM were obtained, guaranteeing the existence of a local 
diffeomorphism from M which identifies V with some partial prolongation of Q . 
However, for applications to integrable systems, nonlinear control theory and 
related areas, it is desirable not only to settle the recognition problem, but ad- 
ditionally, to find a method for explicitly constructing an equivalence between a 
differential system and a partial prolongation, whenever one is known to exist. 
The main aim of this paper is to describe an efficient algorithm for this very 
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construction problem. 

We identify canonical and algorithmically determined integrable Pfaffian systems 
over the ambient manifold whose first integrals determine local equivalences. The 
number of these first integrals that must be constructed, according to our algorithm, 
satisfies a natural minimality criterion. In the special case when V arises from a 
time-invariant nonlinear control system, the complexity of our algorithm agrees with 
that of the well-known GS algorithm [3]. However, even here, the method described 
in this paper enjoys a number of advantages over the former. This point is briefly 
discussed in section 4. We go on in section 5 to give an illustrative example of a 
nonlinear control system where, in addition, a new necessary condition for static 
feedback linearisation is derived. 

Throughout the present work we will rely on the results of the aforementioned 
paper [6] and, wherever possible, adhere closely to its notational conventions. In 
sections 2 and 3, we will recall briefly the salient details of [6] that are required 
to establish the main result which is given in section 4. Section 5 is devoted to 
illustrative examples and an application. 



2. Preliminaries 

2.1. The derived flag. The structure tensor of a sub-bundle V C TM over 
manifold M is a map 

(5o : V X V -> TM/V, 

defined by 

5o(x,y) = po([x,y]) 

where 

Po : TM -^ TM/V 

is the natural projection. If 5o has constant rank on M we define the first derived 
bundle V'"^-* of V as the maximal sub-bundle of TM such that 

Po(V(^)) = Im5o. 

In the constant rank case, we can iterate this procedure, obtaining the sequence of 
structure tensors 

6i: V(') X V^^) ^TM/V^'\ i> 1 

each assumed to have constant rank on M. Define the sub-bundle V^*-* C TM, i > 1, 
as the maximal sub-bundle of TM satisfying pi-iiV^"^') = Im 6i-i where Pi_i : TM 
-^ TM/V^^"^) is the natural projection. The sub-bundle V^*) C TM is the i*^ 
derived bundle of V. We call a bundle V C TM regular if the ranks of all structure 
tensors are constant on M. In this regular case, for dimension reasons there is an 



integer fe > such that V^'^"'"^^ = V^'^^ The smahest such integer is cahed the 
derived length of V and we evidently have a flag of sub-bundles 

Vc V^i) c V(2) c ••• c V^'^) 

called the derived flag of V. Plainly, the numbers dimV'^*^ are diffeoniorphisni 
invariants of V. 

2.2. Cauchy bundles. The Cauchy system or characteristic system x(^) of ^ is 
defined by 

x(v) = {X G r(v) I [x,r(v)]cr(v)}. 

Even if V is regular, its Cauchy system need not have constant rank on M. However, 
if it does, there is a sub-bundle Char V Q V, the Cauchy bundle, whose space of 
smooth sections is xC^^)- 

Definition 2.1. If a regular bundle V C TM satisfies the condition that x(V*^*^) 
has constant rank on M for each z > then V is said to be totally regular. For each 
i, the bundle Char V*-*-* is called the Cauchy bundle or characteristic bundle of V^*-*. 
The elements of Char V*-*^ are called Cauchy vectors. 

It is easy to show that any Cauchy bundle is integrable. 

Definition 2.2. Let V C TM be a totally regular bundle. By the derived type of 
V we shall mean the list of sub-bundles {V*^'-', Char V*^'-*} for all / > 0. 

2.3. The singular variety of a sub-bundle. The final tool we require is less well 
known than the derived type of a bundle, at least in the form in which we shall use 
it, so we discuss this in more detail in the remainder of this section. 

Let V C TM be such that its structure tensor 5o has constant rank on M. 
For the remainder of this section we will not need to retain the subscript 0. For 
simplicity of notation, we will henceforth cease to distinguish between a bundle and 
its module of smooth sections. By PV we denote the projectivisation of V. 

Fix a basis Xi, . . . , X„ for V. If Zi, . . . , Zg complete Xi, . . . , Xn to a frame on 
M, then there are functions c^^ on M, antisymmetric in lower indices, such that 

6{X^,Xp) = cipZk, 1 < a,/3 < n, 1 < A; < s, 

where Z^ is the coset with representative Z^. The elements of PV are lines deter- 
mined by elements e^X^ G V and will be denoted by the symbol E = [e°'Xa\. If E 
is a fixed direction in V then we may seek directions F = [/°Xq] such that 

5{e''X^,fXp)=Q. (2.1) 

Equation (2.1) can be expressed in the form 






^ (e°r-e^r)cS^ = 0, l<k<s, (2.2) 

Q</3=2 
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viewed as a homogeneous linear system for F and in the alternative form 

a{E) • f = 0, (2.3) 

where f = {f^^, f^, ■ ■ ■ , /")• Letting e = (e^, e^, . . . , e") we have clearly that f = e 
is always a solution of (2.3). 

Definition 2.3. The matrix a{E) in (2.3) determined by equations (2.2) will 
be called the polar matrix of the point E G PV. Such a line will be called singular 
if its polar matrix has less than generic rank. The set of all singular lines in V will 
be denoted by the symbol Sing(V). We will denote the set of singular lines in V 
over a point x £ M, by Sing(V)(x). 

Lemma 2.1. For each x € M, Sing{V){x) C PV^, is a linear determinantal 
variety. 

Proof. On the one hand, for each x G M, the polar matrix (j{Ex) has entries which 
are linear functions of the affine coordinates in PV^. On the other, Sing(V)(a;) is 
determined by equating the minors of cr{Ex) to zero. D 

Definition 2.4. The set Sing(V) of all singular points of PV, wih be called the 
singular variety of V. 

We now describe a very useful invariant associated to each point of PV. 
The structure tensor (5 of V C TM induces a map 

degv : PV ^ N 

well defined by 

degv([^]) = dim Image A^, [X] € PV 

where for each X £ V, 



is defined by 



Ax-.V ^ TM/V. 



AxiY)=5iX,Y), YeV. 



Definition 2.5. Let V C TM be a sub-bundle and [X] G PV. The integer 
degv([^]) will be called the degree of [X]. 

Lemma 2.2. For any E G PV, degy{E) is a diffeomorphism invariant. 

Note that if degv([^]) = 0, then X is a Cauchy vector and hence Lemma 2.2 
implies, as a special case, the elementary fact that Char V is an invariant sub-bundle 
ofV. 

Lemma 2.3. For any point E G PV, degv{E) = rank a{E). 
Proof. In fixed bases for V and TM/V, the polar matrix cr{E) is the matrix of the 
vector bundle morphism Ax- D 



Remark 2.1. Amplifying Lemma 2.2, Sing(V) is a diffeomorphism invariant 
in the sense that if Vi,V2 are sub-bundles over Mi,M2, respectively and there is 
a diffeomorphism cf): Mi -^ M2 that identifies them, (/)*Vi = V2, then Sing(V2) 
and Sing(0*Vi) are equivalent as projective subvarieties of PV2. That is, for each 
X € Ml, there is an element of the projective linear group PGL{V2\ ,^^1^) that 
identifies Sing(V2)(0(x)) and Sing(0*Vi)((/)(x)). 

We hasten to point out that the computation of the singular variety for any 
given sub-bundle V C TM is algorithmic. That is, it involves only differentiation 
and commutative algebra operations. In practice, one computes the determinantal 
variety of the generic polar matrix cr{E) as a sub-variety of PV. 

2.4- The singular variety in positive degree. One frequently encounters bundles 
possessing non-trivial Cauchy bundles and it is then natural and prudent to consider 
the notion of degree of a line in the quotient by Char V. Continue with sub-bundle 
V C TM and let 

7?: TM -^ TM/ Char V := TM 

be the natural projection assigning an element of TM to its coset. The structure 
tensor S oi V drops to the quotient bundle V := V/Char V, leading to a tensor 

5: V X V -> TM/V 

well defined by 

S{X,Y) = n{[X,Y]) modV, yX,YeV 

where X = 7r(X), Y = 7r{Y). By analogy, we call 6 the structure tensor of V. As 
before, we are able to introduce a map 

deg9 : PV ^ N 
deg;p([X]) = dim Image A^^ 

A^{Y) = 6{X,Y), V feV. 



well defined by 



where for X £V, 



Here, by [X] we denote the distribution spanned by X € V and Char V. Again, by 
analogy we call deg;p([X]) the degree of [X] G PV. 

Remark 2.2. All definitions and results of subsection 2.3 hold mutatis mutandis 
when the structure tensor 5 is replaced by 6. In particular, we have notions of polar 
matrix and singular variety, as before. Note however, that each point of PV has 
degree one or more. 
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2.5. The resolvent bundle. In [6], we introduced the notion of a Weber structure 
which we now recah. Suppose a totally regular sub-bundle V C TM of rank c+q+1, 
q > 2,c > is defined on manifold M, dimM = c + 2g + 1. Suppose further that 
V satisfies the following additional properties: 

i) dim Char V = c, V^^) = TM, 
ii) S := Sing(V) =¥B^ MP^-i for some rank q sub-bundle BcV. 

Definition 2.6. We will call {V,FB) (or (V, S)) satisfying the above conditions 
a Weber structure on M. 

Given a Weber structure {V,FB), let 7?,g(V) C V, denote the largest sub-bundle 
such that 

n{n^{V))=B. (2.8) 

Definition 2.7. The rank q+c bundle 7^g(V) will be called the resolvent bundle 
associated to the Weber structure (V, S). The bundle B determined by the singular 
variety of V will be called the singular sub-bundle of the Weber structure. A Weber 
structure will be said to be integrable if its resolvent bundle is integrable. 

Remark 2.3. Note that an integrable Weber structure descends to the quotient 
of M by the leaves of Char V to be the contact bundle on J^(]R,]R'^). The term 
'Weber structure' honours Eduard von Weber (1870 - 1934) who, as far as I can 
tell, was the first to publish a proof of the Goursat normal form [7] . 

We record the following properties of the resolvent bundle of a Weber structure. 

Proposition 2.7. ([6]) Let (V, S) be a Weber structure on M and B its singular 
sub-bundle. If q >3, then the following are equivalent 

1) Its resolvent bundle 7?.g(V) C V is integrable. 

2) Each point ofT, = Singiy) has degree one 

3) The structure tensor 6ofV vanishes on B 

Proposition 2.8. ([6]) Let (V, S) be an integrable Weber structure on M . Then 
its resolvent bundle 7^g(V) is the unique, maximal, integrable sub-bundle ofV. 

Remark 2.4. Note that checking the integrability of the resolvent bundle is 
algorithmic. One computes the singular variety Sing(V) = P-B. In turn, the singular 
bundle B algorithmically determines 7^g(V). 

We conclude this section by mentioning some notation. Firstly, in this paper 
we work exclusively in the smooth (C°°) category and all objects and maps will 
be assumed to be smooth without further notice. Secondly, we will often denote 
sub-bundles V C TM by a list of vector fields X,Y,Z, . . . on M enclosed by braces, 

V = {X,Y,Z,...}. 

This will always denote the bundle V whose space of sections is the C°°(M)-module 
generated by vector fields X,Y, Z, . . . . 



3. Partial prolongations and Goursat bundles 

The contact distribution on the 1** order jet bundle of maps from TO — > IR"? 
^(M,M''), q > 1 will 
contact coordinates as 



J"^(M,M'^), q > 1 will be denoted by the symbol Q and locally expressed in 



4'^ = {^- + E^i^-^'^.f}- (3.1) 

A partial prolongation of Q ' may be expressed in contact coordinates as a distri- 
bution on J''{R, W) of the form 

c(r) = {d. + izE £ ^/:;V.--,9,.„,.}. (3.2) 

a = lja=l la=0 

Here and elsewhere in this paper, the symbol r denotes the type of the partial 
prolongation which may be specified by an ordered list of ordered pairs 

T = {qi,ki;q2,k2;... ■,qt,kt) (3.3) 

indicating that there are qa variables of order ka and we use the convention 1 < 
ki < k2 < • ■ ■ < kt] the qa are any positive integers. The positive integer t will be 
called the class of C{t). 

It will often be much more convenient to express the type of a partial prolonga- 
tion as an ordered list oi k := kt non- negative integers 

T= {pi,P2,--- ,Pk) (3.4) 

where the j^^ element pj indicates the number of variables of order j. In this 
notation, if pj = for all j in the range 1 < j < A; — 1, then the contact distribution 
(3.2) has class 1; its type is (0,0, ... ,0,q), q > 1. Such a contact system is a total 
prolongation of Cq , denoted, Cq , an instance of a partial prolongation; namely, 
the (standard) contact distribution on J^(M, M'^). Note that, in the general case, 
the derived length of a partial prolongation of type t is k = kf. 

For any totally regular sub-bundle V C TM, we have the notion of its derived 
type. In section 2, we defined the derived type of a bundle as the list of all derived 
bundles together with their corresponding Cauchy bundles. We shall frequently 
abuse notation by using the term 'derived type of V for the ordered list of ordered 
pairs of the form 

where nij denotes the rank of the j*^ derived bundle V"^ of V and X"' denotes the 
rank of its Cauchy bundle. Char V^-'-'. 
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It is important to relate the type of a partial prolongation to its derived type. 
For this it's convenient to introduce the notions of velocity, acceleration and deccel- 
eration of a sub-bundle. 

Definition 3.1. Let V C TM be a totally regular sub-bundle of derived length 
k and derived type 

'^v = [[mo,x°],K,X^],--- Amk,X% (3.5) 

The velocity of V is the ordered list of k integers 

vel(V) = (Ai,A2,... ,Afc), 

where, 

Aj = rrij — rrij-i, i < j < k. 

The acceleration of V is the ordered list of k integers 

accel(V) = (A2,A2,... ,A|,Afc), 

where 

A- = Ai - Ai_i, 2<i<k. 

The decceleration of V is the ordered list of k integers 

deccel(V) = (-A2, -A^, . . . , -A|, A^). 

(k) 

Note that total prolongations Q ^ have deccelerations of the form 

(0,0,... ,0,g), g>l, 

where there are k — 1 zeros before the last entry q. The classical Goursat normal 
form is the case q = 1 in this family of deccelerations. The main result of [6] is a 
generalisation of this classical result to arbitrary deccelerations. The main aim of 
this paper is to describe an efficient algorithm for constructing the corresponding 
contact coordinates in this general case. 

To recognise when a given sub-bundle has or has not the derived type of a partial 
prolongation (3.2) we introduce one further canonically associated sub-bundle that 
plays a crucial role. 

If V has derived length k we let Char V-i\ denote the intersections V*^-'""'^-' fl 
Char V*--'-', 1 < j < A; — 1. It is easy to see that in every partial prolongation these 
sub-bundles are non-trivial and integrable. 



Proposition 3.1. ([6]) Let sub-bundle V C TM be totally regular with velocity 
and acceleration (Ai, . . . , A^) and (A2, . . . , A|, A^), respectively. Then V has the 
derived type of a partial prolongation C{t) of some type r if and only if 



(3.6) 



mo = l + P, nil = 1 + 2P, 

I 
mi = l + {l + l)P + J2i^ + 1 - J)^P 2<l<k, 

X'' = 2m j — rrij^i — 1, < j < k — 1, 
Xl-i = mi-i - 1, 1 <i<k-l. 

where, P = X]i=i Pi> ^^d 

mi = dimV*^'-', x' = dim Char V*^'-*, Xz-i = dim Char v/i\. 

The type r in C{t) is given by the decceleration, t = deccel(V). 

Definition 3.2. A totahy regular sub-bundle V C TM of derived length k wiU 
be called a Goursat bundle of type r if 

[i] V has the derived type of a partial prolongation whose type is r = deccel(V) 
[ii] Each intersection Char Vl_i is an integrable sub-bundle whose rank, assumed to 

be constant on M, agrees with the corresponding rank in C{t) 
[iii] In case A^ > 1, then V''^"^-' determines an integrable Weber structure on M, 

with singular sub-bundle P^« MP^''-^ 

4. Efficient construction of contact coordinates 

In this section we give our main result, an efficient algorithm for the construction 
of contact coordinates for any smooth distribution locally equivalent to a partial 
prolongation of the contact distribution Q on J^(M,M'^), q>l. Specifically, we 
show how to algorithmically construct certain canonical Pfaffian systems over the 
ambient manifold and appropriate first integrals of these Pfaffian systems whose 
derivatives generate a local equivalence, identifying a given Goursat bundle with 
some partial prolongation. The type of V is given by the decceleration vector, 
deccel(V). 

The main result upon which our algorithm is based is the generalised Goursat 
normal form established in [6]. 

Theorem 4.1. [Generalised Goursat Normal Form], [6]. Let V C TM be a 
Goursat bundle over manifold M , of derived length k > \ and type r = deccel( V) . 
Then there is an open, dense subset M Q M such that the restriction ofV to M is 
locally equivalent to C{t). Conversely any partial prolongation of Cq is a Goursat 
bundle. 
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This theorem settles the recognition problem for a distribution in terms of simple 
constraints on its derived type. While the proof in [6] is constructive, it is extrava- 
gant with respect to the number of integrations that are carried out. The aim here 
is to show that the number of integrations that must be carried out to actually 
construct an equivalence for any Goursat bundle is comparatively small. In fact 
we show that the number of first integrals that must be computed in the general 
case is equal to the number of "dependent variables" featured in C{t), plus one, 
where r = (pi,P2, • • • , Pk) is the bundle's decceleration vector deccel(V). That is, 
^ • -^ pj + 1 = P + 1 first integrals must be found in the general case. The remain- 
ing coordinates are computed by differentiation. This is the "natural" minimality 
criterion alluded to in the Introduction. Our algorithm therefore performs as well 
as the GS algorithm [3] in the special case when distribution V happens to arise 
from a time-invariant control system. 

However, our approach is not restricted to time-invariant control systems, nor, 
indeed to general control systems. Furthermore, it doesn't involve the construction 
of structure equations before an equivalence can be found. Moreover, the geometric 
data that must be computed to settle the recognition problem is expressed more 
naturally in terms of the bundle's derived type. We now procede to the description 
of our method. 

Let V C TM be a Goursat bundle over M of derived length k. Recall that there 
is a distinction between the cases pk > ^ and p^ = 1. In the former case we have 
the filtration 



Char Vi^^ C Char V^^^ C • • • C Char VJi\ C Char V^^^ C • • • 



(4.1) 

k—1 — Sfc-i 



• • • C Char vfZ-^^ C Char V^^"^) C U^ iy^^-^^) c TM 



where T^g (V'^ ^) is the resolvent bundle of the integrable Weber structure 

(VC^-i), Sfc_i). Note that for j in the range 1 < j < /c - 1, Char VJi\ = Char V^-^^ 
if and only if A^_,_^ = 0. 



We will use the convenient notation 



Vj = Aj, 1 < j <k, 

m = iyi+i, 0<i<k-l, (4.2) 

Ni = dimM -mi, <l <k. 



Also, let 



S^^^ = Char V^^^^ -'-' - '-'-'^ 
and 



H(^) = Char V^^^ , e'-^I^ = Char VJi\ , l<j<k-l 



T^,Jv^'-'^) = n^^jv^'^-'^)\ 
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Then we have a filtration of the cotangent bundle T* M 



Tg^_^(V('=-i))cH(^-i) 



c 



:(fc-i) 



C 



c 



:(i) 



C S 



(1) 



CT*M 



(4.3) 



by integrable sub-bundles. It follows easily from Proposition 3.1 and (4.2) that 



diniTg^ fV('=-i)) = iVfe_i + 1, 



diniH^-^') = Nj + A,-+i + 1, dimi 



:(J) 



AT, 



j-i 



+ 1, l<j<k-l, 



and 



dimSj!^! - dimS(^) = -A^.^^ = pj, l<j<k-l 



(4.4) 



(4.5) 



We can therefore construct a filtered basis for sub-bundle H), C T*M as follows 



"(1) 



6Jo,Wl,... ,WAr^._i for Ey}^, 

uJo,uJi,... ,ujNj+nj for E'^^\ I <j <k-l. 

L^o,c^i,... ,cuw._, for Tg^_^(V(^-i)), 



(4.6) 



Definition 4.1. For each j G {1, . . . , A; — 1} such that pj > 0, we refer to the 

sub-bundle Qj = {a;Ar^_|_„^_|_i, . . . ,u!Nj+u-j} C Sj_i as the fundamental bundle of 
order j. 

For each j, for which pj > 0, Qj is a sub-bundle of the integrable bundle 



-^j-i 



and we note the decomposition H^_^ = H^-'^ ©ri^. By the Frobenius theorem, there 
are functions {v''^'''}f^=i and multiplier matrices Mj- with detMj- 7^ such that 

d^^ 
dip^ 



fdip^^'.\ f^N,+n,+l^ 



= M, 



\dipp^'^ ) 



^Nj +nj +2 



\ u:n,+i., / 



mo 



ds(^-). 



(4.7) 



Definition 4.2. We refer to the functions v?'^'-', j G {1,2,... ,k — 1}, defined 
by (4.7) as fundamental functions of order j. Let ^^'^ , ^'^''^, ■ ■ ■ , (p^''-^'^ span the 
first integrals of the integrable sub-bundle Tg (V'*^"^-*). We refer to these as 
fundamental functions of order k. 

We now prove that the construction of the fundamental functions of all orders 
is the only integration that need be carried out in order to construct contact cor- 
rdinates for any Goursat bundle V. 

Theorem 4.2. Let V C TM be a Goursat bundle of derived length k with deccel- 



eration vector t = deccel(V) = {pi, p2, ■ ■ ■ , Pk), Pk > 2. Let {x, p 



l,fc 2,fc 



-v? 



Pk,k 



}| 
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denote the fundamental functions of order k and for each j in the range 1 < j < k — 1 
for which pj > 0, let {^p^'^ , ■ ■ ■ , ip^^'^} denote the fundamental functions of order j 
defined on some open subset U Q M . 

Then there is an open, dense subset U QU and a section YofV such that on lA, 
Yx 7^ and the fundamental functions x, ip^''' := ip^^'^ , together with the functions 

^l:U = ^^, j€{l,...,k}, l<lj<Pj, 0<s,<j-l (4.8) 

are contact coordinates for V, identifying it with the partial prolongation C{t). 

Proof. Fix a point yo & M in a neighbourhood U of which V is a Goursat 
bundle. The proof of Theorem 4.1 in [6] shows that we may extend the fundamental 
functions of order j G {!,... ,k}, for which pj > 0, namely, x = ip^''^, z^^'^ = i^'^'-^, 
1 < ^j < Pj, to a system of contact coordinates 



k,pj 



x,z'-^4-^4■'^... ,4-^');';_ ^ (4.9) 



on an open set lA Q U. Since Char Vq has codimension 1 in V, it follows that 
there is a section Z of V such that Zx = 1 on a dense open subset of ^, which need 
not, in fact, contain y^, and which, for simplicity, we denote by the same symbol. 
Let 

be the local diffeomorphism defined by the change of variable from the original 
coordinates y to the contact coordinates (4.9). Then we have that ^^.V = C{t), 
where 

c{T) = {x = d^+ Y. EE4+i^.'--^.'-^>' 

je{i,...,fc}/,=ih,=o ''^ 

and Zq''' := z^^'K Consequently, we have for some functions a^a^^'^ on hi 



^P^Z = aX+J2o^^"^d^i^ 



U=i 



Since 
we have 



a = {'4}^Z)x = Z{x o V) = Zip^^'^ = 1 

Pj 

V',Z = X + J^a'^'^S^i,,.. (4.10) 
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h >3 ,^h^J ,Jj^^ 



Define functions cp-^' ,{p2' , . . . ,(Pj" for j G {1,2, ... ,k} such that pj > by 

X = ^''\vX:U = z^'j;', i<ij<Pj,o< s, < j - 1, (4.11) 

where ip^''' := ip^^'K By (4.10) we have 

^'{'^ = z(<^'-^-)(y) = zirz'^^^m = iAZKz'^'nwy)) = riz[^''). (4.12) 

The calculation in (4.12) can be repeated so that for each rj in the range 1 < 
fj ^ i, we have (frj'' = '4^*{zrf) showing that the functions defined in (4.11) are 
independent on lA. By their definition, functions (4.11), together with 93*^'^, (/j'^'-^, 
are contact coordinates there. The proof is completed by observing that if Y is 
any section of V such that Yx =/^ onU then we may take Z to be the vector field 
{Yx)-^Y. a 

The only case remaining is pk = 1. Here (M/Char V*-'^~^\ V/Char V*-*^"^-*) is 
a 3-dimensional contact manifold so there is no canonical maximal integrable sub- 
bundle of V^'^"^-' as there is in the case pk > 2. The role of the resolvent bundle 
when Pk = 1 is played by a locally defined bundle, II'^, whose construction we now 
describe. Let x denote any first integral of Char V^*^"^^ and seek any section Z of 
V such that Zx = 1. Define a sub-bundle II'^ C V*^*^"^-' inductively by 

tf+i = [U\Z], n° = Char V^^^ < / < A; - 1. (4.13) 

The proof of Theorem 4.1 shows that II'^ is integrable, has codimension 2 in TM 
and first integral x. In this case, filtration (4.3) is replaced by 

n^^ C H(^-i) C Si^-'^ C • • • C H(i) C S^') c T*M (4.14) 

and the filtered basis (4.6) for Hq is replaced by 

uJo,uJi,... ,ujn,_j for E^^}^, 

ujo,uji,... ,ujNj+nj for E^^\ I < j <k-l, (4.15) 

uo,uJi,... ,uJNk-i for n^ . 

Then by an argument similar to that of Theorem 4.2, we have 

Theorem 4.3. Let V C TM be a Goursat bundle of derived length k and 
decceleration vector t = {pi,p2,... , Pk) , Pk = 1- Let 11^ be the bundle locally 
defined in (4.13). Let 93"*^'^ be any other first integral of\l^ such that dxAdip'^''^ 7^ 
on an open set U C M . 
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Then there is an open, dense subset lA ^lA upon which is defined a section Z of 
V satisfying Zx = 1 such that the fundamental functions x,c/9q""' := (/?'•"■' together 
with the functions 

V^^'li = X^'if, j G {1, . . . , A;}, < Sj < j - 1, l<lj< Pj. 

are contact coordinates for V on U, indentifying it with the partial prolongation 
C{t). 

Theorems 4.2 and 4.3 prove the correctness of algorithm Contact. 
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Algorithm Contact A 



INPUT: Goursat bundle V C TM of derivedlength k and type t = deccel(V) = (pi, ... ,pk), 

Pk > 1- 

1. Build filtration (4.1) of TM and filtration (4.3) of T*M. 

2. Build filtered basis (4.6) of bundle Hq constructed in step 1. 

3. For each j, if pj > 0, compute the fundamental bundle Q,j of order j. 

4. For each j, if f2j is not empty, compute the fundamental functions {(/p j '-'}; :'_-^ of order j. 

5. Fix any fundamental function of order fc, denoted x and any section Z oiV such that Zx = 1. 

6. For each j, if pj > let zh'i = i/^'j'^, 1 < ^j — Pj- 

7. For each j, if pj > define functions 

x,z^'' := zh'J = <^'..J-, :^^^l^ = Zzjf, Q < Sj < j - 1, 1 < I, < p,. 
OUTPUT: Contact coordinates for V identifying it with C{t). 

Algorithm Contact B 



INPUT: Goursat bundle V C TM of derived length k and type t = deccel(V) = (pi, ... ,pk), 

Pk = 1- 

1. Compute filtration (4.1) of TM up to Char V'~^-^\ 

2. Fix any first integral of Char V' -', denoted x, and any section Z of V such that Zx = 1 . 

3. Build distribution 11 , defined by (4.13), giving refinement (4.14) of the filtration constructed 
in step 1. 

4. Let z^ := ^p^'^ be any first integral of H*^ such that dx A di^^^^ /; 0. 

5. Build filtration (4.14) of T* M from that of TM constructed in steps 1 and 3. 

6. Build filtered basis (4.15) of bundle Sq . 

7. For each j, if pj > 0, compute the fundamental bundle flj of order j. 

8. For each j, if Q.j is not empty, compute the fundamental functions {(/? J '•'}( ^_-^ of order j. 

9. For each j, if pj > let z'^i'^ = v^'^'^, 1 <^ ^j ^ Pj- 
10. For each j, if pj > define functions 

x,z^'' := zh'^i = <^'..^ 2^^;:i^ = Zz^^\ < Sj < j - 1, 1 < /j < Pj. 

OUTPUT: Contact coordinates for V identifying it with C{t). 



5. Examples and Applications 

In this section we illustrate our algorithm for finding contact coordinates for dis- 
tributions locally equivalent to partial prolongations and indicate some applications 
to nonlinear control theory. 

Example 5.1. We begin with a simple illustrative example that is sufficiently 
non-trivial to illustrate the all the steps in algorithm Contact A. Consider the 
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sub-bundle of TM^^ defined by 

V = le'^^^dx^ + (a;4 - xs - e'^^)dx2 + (^5 - 2:21 - 2:20 + X4 - ais - 6^1)9^3 + 
(x9 + 2x2i +22:20 - 4x4+ 42:8 + 3e^i) (9^4 + d^g) + (2 - 2xi2 + XT)dxe + 

{xs + e^^)dxj - dxg, + (9x9 + (1 - 2:13 - Xl4 - Xl2)(9xio + 

(-(a;i3 +a;i4) + X4 - xg - e^i)^^;!! + -(a^is - 2;i4)92;i2 + 

(xi5 - xi6)9xi3 + (a;i5 + xif,)dx^i + xitS^^j^ + (xis - 2x2)dx^e + 

(xi9 - X8 - e^i )93;i7 + 3;2o9xi8 + (3x21 + 8x20 - 6x4 + 6x3 + 5e^i + X2)dx^g , 

dx4 + 29x20 1 ^^5, 9x9, ^^20 ~ '^a;2n C^X5 — 29x9 + '^2:21 \ = {-'^1,^2, • • • i-'^e}- 

We compute tlie (refined) derived type to be 

[[6, 0], [f 1, 5, 7], [14, 10, 10], [17, 13, 14], [19, 16, 16], [21, 21]] 

where for i = and i = 5 the two-element lists record [w,i,x*] and for 1 < i < 4 
the 3-element lists record [?7ii, Xi-n X*]- Hence the derived length is fc = 5 and by 
Proposition 3.1, we confirm that V has the derived type of a partial prolongation 
whose type is (2,0,1,0,2). Since ^5 = 2 > 1, we check the singular variety of 
9(4) ^ v^'^VChar V^"^). We compute that 

V(") = {Xi = ^(Xi), X2 = ^(9,,J, X3 = ^(a,,3 + a,,J} 
whose non-zero structure is 

?(Xi,X2) = {9{2d,,+d^J), 6{X^,Xs) = {^{2d^,,-d^,,)). 
Here 5 is the structure tensor of V'^'^-', 

^ : TM^i -^ TR^i/Char V^"^) =: TM^i 

is the natural projection and for Z G TM^^, (tt{Z)'^ denotes the element of the 
quotient bundle 

fM2i/v(4) 

with representative Tr{Z). From this, we easily compute that an arbitrary point 
E = [a^Xi + a^Xs + a^Xg] G PV^"^) has polar matrix 

f -c? a" \ 
V-a3 aV' 

from which we deduce that 

sing(v(^)) = p{5?(a,,j,5?(9,,3 + a,,j} 
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and the resolvent bundle is therefore 

7eg^(V(^)) = Char V^^^ {d^,,,d^,, + 9,, J. 

The filtration induced on T'M?^ by V is 

Char V^^^ C Char V^^^ C Char V^^^ cChar vf ^ C Char V^^^ 

(5 1) 
C Char V(^) C 7e£^(V(4)) C TM^i 

and one can check that the resolvent bundle is integrable. This data confirms that 
V is a Goursat bundle of type deccel(V) = (2,0,1,0,2). By Theorem 4.1, V is 
locally equivalent to the partial prolongation with this type. We now go on to use 
algorithm Contact to construct an explicit equivalence. 

The filtration of T*^"^ induced by V and dual to (5.1) is, in this case, 

Tg^(V(^)) C S(4) c H(3) c Hf C S(2) c H(i) C ~\l^ C T^M^i (5.2) 

Pursuing step 2 of algorithm Contact A, a filtered basis for Hq is 

Tg (V^'^-') = {dxi,dx2 — (ixii,dxio} 

rS ' = {dxi, dx2 — dxii, dxiQ, dxi2, dxi^ + dxi^} 

'E^'^' = {dxi, dx2 — dxii, dxiQ, dxi2, dxi^ + dxi4, dxi^, dxi^} 

^2 = {dxi, dx2 — dxii, dxio, dxi2, dxis + dxi^, dxi^, dxi4, dx^} 

S'^-* = {dxi, dx2 — dxii, dxio, dxi2, dxi^ + dxi4, dxi^, dxi4, dxQ,dxT, dxiQ, dxn} 

S'^-' = {dxi, dx2 — dxii, dxio, dxi2, dxis + dxi4, dxi^, dxi4, dxQ,dxT, dxi^, dxn, 

dxiQ,dxs,2dxii — dxis} 

Sg = {dxi, dx2 — dxii, dxio, dxi2, dxis + dxi4, dxi^, dxi4,dxQ,dxT, dxiQ, dxn, 

dxig,dxs,2dxii — dxis, dx^, dxis} 

Since pi,P3 and p^ alone are non-zero, we deduce from this that the fundamental 
bundles are 

Qi = {dxs,dxis}, ^3 = {dxa}, O5 = {dxi,dx2 — dxii,dxio} 

and the rest are empty. The corresponding fundamental functions may therefore 
be taken to be (for instance) 

•^1 = {z^'^ = X3,z'^'^ =Xi8}, J^3 = {Z^'^ = Xq}, J^5 = {z^'^ = -X2+X11, z"^'^ = Xio}| 
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and we take the independent variable x to be e^^ , for then Xie^^ = 1. Consequently, 
Z = Xi may be taken to be the operator of total differentiation. 

Finally, executing step 7, we compute the remaining components of the equiv- 
alence by differentiating once the elements of jFi, differentiating three times the 
elements of J-^ and finally differentiating five times the elements of J-^ : 

Zi — AiZ , ^^3 + 1 — ^l^ra 1 ^rs + l " ^l^rs 

where 

1 < ^1 < 2, 1 < /s < 2, 

< S3 < 2, < S5 < 4 

and Zq' := z^'^ . We thereby obtain the functions 

X3,z'^'^ = Xis,Zj^' = X5 — X21 — 3:20 + X4, - Xs — e'^^ , Zj^' = X20, 

1.3 li3 n I I ri 1,3 I I X^ 

z' =XQ,z^' = -2x12 +xr + 2,Z2 = 2;8 — xis +2;i4 +e 1, 

Zg' =2:9+ 2(a;20 + 2:21 - 22:4 + 2x8 + ^w) + Se^^i , 
z^'^ = 2:10,2;^'^ = -a;2 + xii,z^' = 1 - (2:12 +X13 +2:14), z^' = 1/2(2:13 +2:14), 
^2' = 1/2(2:14 - 2:13 - 4x15), 2:2' =xi5,Z3' = xi6 - 2x17,^3' =2:17, 
z^' = xi8 - 2(x2 + xi9 + X8 +e^i),2:4' = xig — X8 - e^i , 
z^' = 2(x4 - X8 - e^i - X21) - X2o,25''' = 2(-X4 + X8 + e^i) + X20 +X21, 

in accordance with Contact A. By Theorem 4.2, these functions define a local 
diffeomorphism ip: M?^ -^ j(2,o,i,o,2) satisfying 

V'*V = C(2,0,1,0,2) 

where J\2:0,i,o,2) jgjiotes the partial prolongation of J^(M,M^) in which two vari- 
ables remain at order 1, one variable is prolonged to order 3, two are prolonged to 
order 5. Finally, the contact distribution on J\2:0:i:0,2) j^^g ^-j^q form 

2 2 2 4 



C(2,0,l,0,2) = {d, + Y^ z[^''d,.,. + Y^ zlf^.d^.^.s + Y.I1 'hf+Ar^ 

9^1,1,5^2,1 5^1,3, d 1,5 82, 5}, 

^1 ^1 ^3 ^5 ^5 

using the notational conventions in the proof of Theorem 4.2. 
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Example 5.2. [Nonlinear Control Theory]. In geometric control theory one 
considers a car moving in the xy-plane (see [8]) modehed as foUows. The state 
of the car is described by four variables {x,y,6,(p). The ordered pair {x,y) gives 
the coordinates on the xy-plane of the centre of the rear axle. The variable 9 is 
the angle between the x-axis fixed on the plane and the vertical axis V of the car, 
running perpendicular to the axles. The variable (p is the angle the front wheels 
make relative to V. 

Assuming the wheels do not slip as the car moves in the plane, then one obtains 
the Pfaffian system 

— sin 6 dx + cos 9 dy = 0, L cos (j) d9 — sin (f) (cos 9 dx + sin 9 dy) = 0. 

This in turns leads to the control system / = {oji = 0, . . . , u;4 = 0}, known as the 
kinematic car, where 

ooi = dx — v^ cos 9 dt, ll>2 = dy — u^ sin 9 dt, cos = d9 — tan (f) dt, UJ4 = dc/) — u^ dt, 

ij 

and L is the length from the rear to the front axle. Here u^ models the speed of the 
point (x,y) and v? the speed at which the front wheels swivel. Variables u^^u^ are 
the controls, for prescribing these as functions of time t gives a system of ordinary 
differential equations for the state of the car. 

Equivalently, control system / defines the sub-bundle /C = J"*" C T (Mj x M), 
given by 

^ = {dt + u^ (cos 9 dx + sin 9 dy + — tan (f) dg) + u^ dfp, du^ , du-2 } 

where M = M3 -, xM.J -, x S^ x S^ is the manifold of states and controls. An 

(x,y) (ui,U2) 

important question in nonlinear control theory is: when can a nonlinear control 
system be "linearised" by a static feedback trans forraation and more generally a 
dynamic feedback transformation! We will study these questions for the example 
of the kinematic car as an application of our general algorithm Contact. 

We begin by showing that K, is, in fact, a Goursat bundle of type r = (1,0, 1). 
The refined derived type of /C is 

[[3,0], [5, 2, 3], [6, 4, 4], [7, 7]]. 

Consequently the derived length is 3 and deccel(/C) = (1,0,1). The filtration in- 
duced by /C is 

Char /C|)^) C Char K^^^ C Char /C^^) c TM^ 

where 

Char ICq^ = {9„j,9„J, Char /C*^^^ = {5„,, 9„2, 9t}, Char /C^^^ = {9„j,5„2, ^t,^^}. 

(5.4) 
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Since the distributions in (5.4) are all integrable and p^ = 1, we conclude that /C is 
a Goursat bundle of type (1,0, 1). By Theorem 4.1, /C is locally equivalent to the 
contact distribution C(1,0, 1). This settles the recognition question for /C. We go 
on to find an equivalence using Contact. 

For this, according to Contact B, we compute distribution XI"^. The invariants 
of Char JC^"^' are x, y, 6. Any one of these may be taken to be the independent 
variable. If we choose x for this purpose then we take 

X = dx ~\ — -, T:(dt + u^ sin 6* dy + L~^u^ tanc/) de + u^ ds) 

u^ cos 

for the total differential operator. By (4.13), we find that 

n^ = {du,,du2,dt,d^,de}. 
The filtration of T*M^ induced by /C is therefore 

n3^ c s(2) c s(i) c Ei'^ c r*M7 

where 

n^^ = {dx,dy}, S(2) = {dx,dy,d9}, S^^) = {dx , dy , d9 , dcf)} , s{,^^ = {dx,dy,(ie,#,dt}.| 

It follows that the fundamental bundles are f^i = {dt} and f^a = {dy} and the 
fundamental functions are therefore J-'i = {t}, J-3 = {y}. Consequently, by Contact 
B the map ^: M^ ^ j(i:0,i> defined by 

X, z''' = t, zl^^ = Xt, z''^ = y, z'/ = Xy, z"/ = X^y, z"/ = X^y (5.5) 

is an equivalence, according to Theorem 4.3. That is, ^^./C = C(l, 0, 1). 

While (5.5) is certainly an identification with a partial prolongation, it is not of 
much use in control theory. This is because it is not a feedback equivalence. That 
is, (5.5) does not respect the special role played by the time coordinate t in the 
original control system, nor the distinction between the roles played by the state 
and control variables. We pause briefly to describe the class of transformations that 
preserves the set of all control systems. 

Let 

^=/(t,x,u) (5.6) 

be a control system, where x denotes the state variables and u the control variables. 
A local diffeomorphism of the form 

t^t, X 1-^ $(x), u t^ ^(x, u) (5.7) 
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is said to be a (static) feedback transformation. A question of interest is: does there 
exist a static feedback transformation that identifies a given control system with 
some partial prolongation of Q for some q or, as it is more commonly known in 
the control theory literature, a Brunovsky normal form? 

Writing, (5.5) out explicitly shows that it is not a static feedback equivalence. 
In fact, an elegant result of Sluis [5], applied to the kinematic car system /, proves 
that no feedback equivalence exists. Sluis' theorem relies on the GS algorithm and 
consequently applies to time-invariant control systems such as the kinematic car. In 
fact, a stronger necessary condition for static feedback equivalence can be derived 
via Theorems 4.2 and 4.3. 

Theorem 5.1. Let K. = {dt + /(t,x, u)9x,5u} C T (Mt x M) arise from any 
smooth control system (5.6), where M is the manifold of states and controls. A 
necessary condition in order that /C he static feedback equivalent to a Brunovsky 
normal form is that 
[i] K, is a Goursat bundle 

[ii] // k is the derived length of K, then dt € Char ]C^^~^^ if Pk = 1 or dt € 

r^,JJC^''-'^)tfPk>i. 

Proof. Suppose there is static feedback equivalence i? identifying /C with some 
Brunovsky normal form. Every such normal form is a partial prolongation C{t) 

of Cq , where t = {pi, p2, . . . , Pk), Pk ^ ^- In the case pk > I, the independent 
variable x of C{t) is an invariant of its resolvent bundle. Consequently, 'Q*x = t 
is an invariant of the resolvent bundle TZY,^_-i^{K,^^~'^') determined by /C. In the 
case Pk = 1, X is an invariant of the Cauchy bundle Char C{t)^^~'^' of the k ^ 1 
derived bundle C{tY^~'^\ Hence, 'd*x = t is an invariant of Char K,^^~^\ D 

While for the kinematic car, /C is certainly a Goursat bundle, it does not satisfy 
condition [ii] of Theorem 5.1. We deduce that K, is not static feedback linearisable, 
in agreement with Sluis' theorem. 

However, as is well known for this example, a certain Cartan prolongation of K, 
is static feedback linearisable. We will not give precise definitions for this, referring 
the interested reader to [5] for the details. We merely want to present a Cartan 
prolongation of K, and show how to apply algorithm Contact to determine a static 
feedback linearisation of the Cartan prolonged distribution. 

We obtain a Cartan prolongation of the kinematic car system as follows. Define 
a new control system 

dx 1 dy 1 ■ a dO u^ dcp 2 d^^ 2 dw^ ,^ on 

— - = n cos t>,—-=u sm t/, -— = — tan (p, —- = u , —— = w , — ; — = Vi (5.»j 
dt dt dt L dt dt dt 

by "twice differentiating u^". In control system (5.8) the coordinate u^ has become 
a state variable and the new control variables are v^ and u'^. To see the significance 
for control theory of this admittedly ad hoc construction one needs to check for the 
existence of a static feedback equivalence for system (5.8). We begin by changing 
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notation slightly and examining the sub-bundle pr /C C T (Mj x M), defined by 
(5.8), where M is the manifold of new states and controls. Setting u^ = w^ and 
u'^ = v'^, we have 

pr /C = {T = dt+w^{cos9 dx+sm9 dy^ — d9)+v'^ d^i^+w'^dwi+v^ 8^2, 8^1,8^2}. 

The refined derived type is 

[[3,0], [5, 2, 2], [7, 4, 4], [9, 9]]. 

Proposition 3.1 shows that this is the derived type of the partial prolongation 
C(0,0,2), that is, the total prolongation Q . A calculation reveals that pr /C*^^) := 
pr JC^^' /Char pr JC^"^' is spanned by 

pr /C*^^-* =Tr{T = 8t + w^^cosO 8^ + sin 9 8y), 8^1,80} 

and that the polar matrix of a line E = [a^7r(r) + a^7r(5^i) + a?Tr{8g)] is 

-a^ a^ 
-a^ a^ 

with respect to the non-zero structure 

?(^(a^i),^(r)) = (^(cos^a^+sin^aj,)), ?(^(r),^(ae)) = (5?(z^(sin^9^-cose9j,))).| 

Here, again, 5 is the structure tensor of pr K,^"^', 

vr : T{^t X M) -^ T(Rt x M)/Char pr /C^^) =: f (Mt x M) 

is the natural projection and for Z € T[Rt x M), (7r(Z)) denotes the element of 
the quotient bundle 

f (Mt X M)/pr £(2) 

with representative Tr{Z). It follows that Sing(pr IC^^^) = ¥{Tr{8w^),7r{8g)} and so, 
the resolvent bundle of pr JC^"^' is 

7^g^(pr/C(2)) = Char /C(2)© {3^1,94, 

which is integrable. The above data shows that pr /C is a Goursat bundle of type 
(0, 0, 2). By Theorem 4.1, there is an equivalence identifying pr /C with the contact 
distribution C(0, 0, 2), that is, a Brunovsky normal form. But is it a static feedback 
equivalence? To find out, we use algorithm Contact. 
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The only non-empty fundamental bundle in this case is the one of order 3 

Tg^(pr /C(2)) = 7^£^(pr K.^^'^t = {dt,dx,dy}. 
So in this case, dt G Tg (pr /C^^-*) and we may choose z'^''^ = x, z^'^ = y and 

1,3 _ rp 1,3 2,3 _ rp 2,3 1,3 _ rp 1,3 2,3 _ rp 2,3 1,3 _ rp 1,3 2,3 _ rp 2,3 

Z-^ — -L Z , Z-^ — 1 Z , ^2 — ^ ^1 ; ^2 — ^ ^1 ' ^3 — -^ ^2 ' % — -^ ^2 • 

It is a simple matter to verify that these are indeed contact coordinates as predicted 
by Theorem 4.2. Moreover, an inspection of the formulas 

1,3 1 a 2,3 1 ■ a 

z^ =w cost/, z^ =w siiity, 

1.3 2 n ("'^)^ ■ d+ A 2,3 2 ■ fl , (^^^ n. i, 

22 =10 cost/ sinytan0, ^2 = '"' siii^H cost/ tan 0, 

13 — (ui^)^ cos sin — 3w^w^Lsin(jf)cos()f)sin0 — {w^)'^v'^Lsind -\- v^L^ cosflcos^ (^ 

""3' " L2 C0S2 </> 

2 3 — (ui^)^ sin 6 sin^ </> + 3iLi^iu^Lsin(/)Cos0cos0 + {w^)^v^L cos 6 + u^L^ sinflcos^ (fi 

L2 cos2 

shows that they define a static feedback equivalence. Because of the projection 
M ^ M, any integral manifold of pr /C maps to a unique integral manifold of 
/C. This, together with the fact that there is a static feedback equivalence for 
pr /C implies that the kinematic car example is dynamic feedback linearisable in the 
language of nonlinear control theory. It is an important and largely open problem to 
geometrically characterise the class of nonlinear control systems which are dynamic 
feedback linearisable. 
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